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distributions are applied. For the lower bound estimates, the characteristic function method is 
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1. Introduction 

The closeness of a compound Poisson (CP) distribution to the Markov binomial (MB) 
distribution has been investigated in numerous papers; see, for example, [5, 6, 12, 14, 20, 
31, 34, 37] and the references therein. Related problems were considered in [4, 7, 8, 13, 
18, 19, 30, 33] and [38]. One would expect the MB-CP case to have been comprehensively 
studied. As it turns out, this is not the case. Many papers deal with the convergence facts 
only. Only a few of the papers dealing with the estimates of accuracy of approximation 
involve no assumptions about the stationarity of the Markov chain. 

The aim of this paper is to discuss some compound approximations for non-stationary 
Markov chains. We show that for our version of the MB distribution, the natural approx- 
imation is a convolution of CP and compound binomial distributions, both having the 
same compounding geometric law. We outline some principles of construction of asymp- 
totic expansions and consider second order approximations. Part of the paper is devoted 
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to signed compound Poisson approximations which can be viewed as the second order 
expansions in the exponent. We obtain upper and lower bound estimates and show that 
under certain conditions, they are of the same order of accuracy. All estimates are proved 
for the total variation, local and Wasserstein norms. For the upper bound estimates, we 
employ a convolution technique which can be dated back to [23]. For the lower bound 
estimates, we use the characteristic function method. The methods of proof do not allow 
for reasonably small absolute constants. However, in special cases, asymptotically sharp 
constants are calculated. 

We now introduce some notation. Let Ik denote the distribution concentrated at an 
integer fceZ, the set of integers, and set I = Iq. In what follows, V and M denote two 
finite signed measures on Z. Products and powers of V and M are understood in the 
convolution sense, that is, VM{A} = J^fcL-oo V{A — k}M{k} for a set ACZ; further, 
M° = I. The total variation norm, the local norm and the Wasserstein norm of M are 
denoted by 



|M||= J2 I^WI' 



IMI 



h— — oo 



= sup|Af{fc}|, 



\M\ 



\v 



= J2 |M{(-oo,fc]}|, 



h—— OO 



respectively. Note that ||(7i — /)M||\y = ||M||. The logarithm and exponential of M are 
given, respectively, by 



00 (— Fi fe + 1 1 
lnM = EMfe ( M " ^ (if ll^-^ll < 1), e M =exp{Af} = ]T-M fc . 



k=l 



fe=0 



Note that 



llVAflloo^HVIIIlAfllao, ||VM||<||V||||M||, ||e M ||<ell M H 



Let M(t) (t £ R) be the Fourier transform of M. We denote by C positive absolute 
constants. stands for any finite signed measure on Z satisfying ||0|| < 1. The values of 
C and can vary from line to line, or even within the same line. Sometimes, to avoid 
possible ambiguity, the C's are supplied with indices. For i£R and = {1, 2,3,.. .}, 

we set 



= ^rx(a;-l)---(a;-A; + l), 
k\ 



= 1. 



Let £o)£i! • • • j£n: • • • be a Markov chain with the initial distribution 

P(6> = l)=Po, P(G> = 0) = l-po, PoG[0,l] 
and transition probabilities 

P(^ = l|6-i = 1) =P, P(6 = 0|6-i = 1) = q, 

p(& = i|6-i = 0) = q, P(6 = o|6-i = 0) =p, 

p + q = q + p=l, p,?e(0,l), tGN. 
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The distribution of S n = £1 -| — • + £ n ( n G N) is called the Markov binomial distribution. 
We denote it by F n . We should note that the definition of the Markov binomial distri- 
bution varies slightly from paper to paper; see [12, 30] and [36]. We choose the definition 
which, on the one hand, contains the binomial distribution as a special case and, on the 
other hand, allows comparison to the Dobrushin's results. Dobrushin [12] assumed that 
Po = 1 and considered S n -i + 1- 

Later, we will need various characteristics of S n . Let 

<N qq 2 ( , q \ 7i 

7i = — — j 72~ 



h i h {q + q) 2 

73 =7i73, 

7i. i f 2 - , pq(Zq-q) , Zqq 2 \ , g 

73 = it + , ; _^ i P 9 H — V , : _ N9 f + — — P - 



3 + g + g (q + q) 2 } q + q\ q + 

A = ?i-po, ^i=7i — — -Po , ^2=Po — — , Ci = In— =0.4567. 



.q + q J q + q 19 

We use the following measures also: 

qe lf 




G = qhJ2P Jl J (G(t) 
^i = (l-7iK + 7iG, i/ 1 A =cxp|A^^l7i(G-^j (H*(t) = (H^t))*) 

D j =exp\Y l 'Vi(G-T) i \, 1<j<3, Z?^ cxp{A 7 i(G - /)}. 



2. Known results 

In this section, we discuss some of the known results on the compound Poisson approxi- 
mations to the MB distribution. Many papers deal with the convergence facts only; see, 
for example, [14, 19, 20, 36]. Usually, the chain is assumed to be stationary. A typical 
example is Theorem 4.1 in [31] which states that if po = q/(q + q) and a > 0, then 

\\F n - exp { a(G - 7)}|| < 2\nq- a\ + gO±£+g^) . (i) 

q + q 

Even if we choose a = nq, the order of accuracy in (1) is not better than nq 2 . A similar 
estimate was obtained in Theorem 5 of [37] . If we use the terminology of the book [2] , 
we can say that the estimate (1) contains no 'magic' factor. If we turn to the papers 
with 'magic' factors, then we have the following results. In [5], it was proven that if 
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0<p<C < 1, then 



IK,. 



<Cmax(p ,g)min 1, -= 



Cmin((7, nq 2 ) + Cc 



-Cn 



(2) 



The accuracy can be improved, by some asymptotic expansions, to 



F n -D? (I+po 



q 2 {p-q) 
{q + q) 2 



h 



< Cq(p + q) min ( 1, — = 
yjnq 



+ Cmm(q,nq 2 ) + Cc 



-Cn 



(3) 



Note that in [5] formulas (4.5), (4.12) and (4.23) contain misprints. The parameter po is 
misplaced and should be in the brackets. If p — >• p = const, q — > and nq — > oo, then the 
order of accuracy in (2) is max(<?, (nq)~~ 1 / 2 ). Also, the order of accuracy in (3) is q. We 
can hardly call (3) the second order expansion since the improvement of the accuracy 
was achieved due to the more precise approximation of the initial distribution of £o only. 

The main idea of signed CP approximations is to leave more than one factorial cumu- 
lant in the exponent. In short, the signed CP measure has the same structure as the CP 
measure, but can have negative Poisson parameters. Such approximations are commonly 
used in insurance models and in limit theorems; see [1, 9, 17, 21, 22, 25, 28] and the 
references therein. For the MB distribution in [5] the following result is proved. If 



p<C<l, 



q 



< 



l -C 
30 ' 



then 



\F n -D%\\<C(p + q)lmm 




max(po, q) min 1 



/nq 



+ e 



-Cn 



(4) 



(->) 



Note that 72 < and, therefore, D2 is a signed measure rather than a distribution. As 
a rule, signed CP approximations are more accurate than CP approximations. Indeed, 
if n — > 00, then (5) gives the estimate converging to zero, even if p and q are constants. 
Also, (2) and (3) are non-trivial even if we only have q = o(l). 

We are unaware of any lower bound estimate for the compound Poisson approximation 
to the Markov binomial distribution. 



3. The main results 
3.1. Geometric expansions 

Before formulating our results, it is necessary to explain the choice of approximating 
measures. Dobrushin [12] proved that if p — > p, nq — > A and po = 1, then the limit dis- 
tribution for 1 + S n -i is the convolution Giexp{A(Gi — /)}, where G\ is a geometric 
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distribution with parameter p, that is, G\{t) = (1 — p)c li /(l — pe 1 '). This suggests that 
approximation of S n for arbitrary pq should also be based on expansions in powers of 
G — I. 

Let F be concentrated on Z and have all moments finite. We can write formally (i.e., 
without investigating conditions needed for the convergence of series) 



= 1 + E % eit - !) J = 1 + E % g w - ir 

z — ' ?! z — ' ml 

j = l ■' m=l 



Here, Vj (j = 1,2,...) are factorial moments of F, and v m can be called the geometric 
factorial moments. Since 



it - ?(G(t)-l) 
e — 1 — 



l+ P (G(t)-l) 
it is not difficult to establish a relation between i/j and v n 



j 

j\ \ p J \m — j 



(-^EtH--) (Zz)), 



Similar relations hold for factorial cumulants and geometric factorial cumulants. For the 
MB distribution, we have 

Vi = qv\ = qES n = nji + xi+x 2 - {x\ + x%){p - q) n ■ (7) 

(For the formula of the mean, see [7].) Since we will assume p and q to be small, the 
last summand in (7) will be neglected. As it turns out, F n is close to some convolution 
WiA"; see (32) below. We use (6) for choosing the approximating measure for W\. The 
cumulant analog of (6) is used for A". 



3.2. Compound Poisson approximation 

In this paper, we usually assume that 

P<h ?<M- (8) 

The size of the absolute constants is determined by the method of proof. We expand 
F n as a series of convolutions of measures. The remainder term is usually estimated 
by a series containing powers of q + p. If the sum q + p is sufficiently small, the series 
converges. Thus, although we have some freedom in the choice of magnitude of p and 
q, the sum p + q must be small. The choice of condition (8) is determined by the fact 
that the CP limit occurs when nq— > A; see, for example, Table 1 in [12]. Therefore, if we 
expect the CP approximation to be accurate, then q should be small. On the other hand, 
we have included the case q = constant, that is, the case which is usually associated with 
the normal approximation. We choose the assumption p < 1/2 instead of (4), in order to 
make our proofs clearer. 
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Theorem 3.1. Let p < 1/2. We then have 

\\F n - HDi\\ < Cq(p + q) minf 1, -jL ) + Cmm(q,nq 2 ) + C(p + q)e- Cin . 



Inq 



If, in addition, q< 1/30, then 



(10) 



\\F n - HD^ <Cq(p + q)mm(l,-^j + Cmin\^J ^,nq 2 J 
+ C(p + q)c- Cin , 

\\F n - HD$\\ W < Cq(p + q) + Cmm(q^l,nq 2 ) + C(p + q)e- Cin . (11) 
Corollary 3.1. If (8) is satisfied and nq>l, then 

\\F n -HD?\\<Cq, WFn-HDfWnKcJl, \\F„ - HD[ l \\ w < Cq^riq. 

V n 

Corollary 3.2. If (8) is satisfied, then 

\\F n -HH^\<C(q+pe- c ^). 

Remark 3.1. 

(i) H is compound Poisson distribution; see Lemma 5.3 below. If q < p, then is 
also a CP distribution. Thus, we see that there exist quite different forms of CP 
approximations with similar orders of accuracy. 

(ii) The estimate (9) is slightly better than (3) for p,q< exp{— Cin}, and more accu- 
rate than (2) for q < l/^/ri^ and p n > C. 

(iii) For the closeness of F n and HiD^, it suffices to assume q— > 0, in considerable 
contrast to riq — >• A, the latter being needed for the convergence to the limit CP 
law. 

(iv) We can write min(g, riq ) = riq 2 min(l, (riq) -1 ). The last factor, in terms of [2], 
page 5, can be called the 'magic' factor. 

The accuracy of approximation can be improved by the second order expansion. 

Theorem 3.2. If p< 1/2, then 

\\F n - HD\{I + n l2 (G - I) 2 )\\ < Clq 2 + pqmm (l, -±= ) +(p + q) C - C ^ n 



'nq 



If, in addition, q< 1/30, then 

\\F n -HD*(I + n l2 (G-I) 2 )\\ 
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< C\q 2 min( 1,— = ) + pgmin( 1, — ] + (p + <7)c~ Cin 
\\F n -HD*(I + n l2 (G-I) 2 )\\ w 



< C{q 2 max(l, y'ng) + M+ (p + g)e- Cl ™}. 

Note that the last estimate contains max(l, y/rCq), reflecting the fact that the estimates 
for the Wasserstein distance are less accurate than the ones for the total variation norm. 
It is even more evident when nq > 1 . 



Corollary 3.3. // (8) is satisfied and nq> 1, then the estimates in Theorem 3.2 

Cq(q + ^\ cMq + 4=), Cq^q(q + P 



are 



\ v n 1 J V n \ v n Q ) \ v nc l 

respectively. 

We see that, in general, even the second order estimates in total variation are only 
meaningful for q = o(l). 

3.3. Signed compound Poisson approximations 

The choice of a signed CP approximation, in general, means that the first term of the 
asymptotic expansion, unlike Theorem 3.2, is in the exponent. 

Theorem 3.3. If condition (8) is satisfied, then 

||F n -ifexp{xi(G-/)}^||<a(p + g)| m in^,y|j + C - Cl "j, (12) 

||F„ - H exp{^ (G - I)}D$\\ x < C(p + q) jmin U ^) + c- Cl "| , 
WFn-HcxpiMG-I^DZWw^Cip + lDiq + c- ^}. 

Note that for nq < 1 and po = constant, (12) is more accurate that (5). More impor- 
tantly, when p = constant and q = constant, the estimate (12) is of order 0(n -1 / 2 ). In 
this sense, the signed CP approximation is comparable to the normal one and, moreover, 
it holds in the total variation metric. Meanwhile, for discrete distributions, the normal 
approximation holds in the uniform metric only. Just as in the CP case, the second order 
expansions can be used. 

Theorem 3.4. // (8) holds, then 



\\F n - Hexp{MG - I)}D%(I + n 73 (G - I) 3 )\ 
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(13) 

< C(p + q) I min («,-) + c~ Cl " 



(14) 



(15) 



||F n - HeMMG - I)}D%(I + nj 3 (G - I) 3 ) 

<C(p + q)\mm( q,—^= ) +c-°' in 
I \ n^/nqj 

\\F n - Hexp{MG - I)}D%(I + n l3 (G - lf)\\ w 
<C(p + 3)jminUy^n +e- c ^. 

Corollary 3.4. Let nq>l. The estimates (13)-(15) are then at least of order 
C(p + q) C(p + q) C(p + q)yM 

respectively. 

In Theorem 3.4, only a part of the asymptotic expansion is in the exponent. Therefore, 
the following question naturally arises. Is it possible to find a signed CP measure which, 
up to a constant, provides the same accuracy as in Theorem 3.4? As it follows from the 
following result, such a measure indeed exists. 

Theorem 3.5. If (8) holds, then 

1 



li^-ifexplxiCG -I)}Dl || < C{p + q)\ min (q, 



n 



e -c in 



|f n -ffexp{xi(G-.0}I?J||oo < C(p + q){mm(q,—^= ) + e" c 



\F n - Hcxp{MG ~ I)}D%\\ W < C(p + g) j min (q, JT) + c- Cl " 



3.4. Lower bound estimates 

In this section, we show that in some cases, the estimates in Theorems 3.1 and 3.3 are 
of the correct order. We concentrate our attention on the case riq> 1. 

Theorem 3.6. Let condition (8) be satisfied and let nq>l. Then, for some absolute 
constants C% and C3, 



\F n -HDl\\>C 2 q[l-C 3 [q+-j=) ), 
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\F n - HD1\\ W > C 2 q^/nq[ l-C 3 [q 



nq / 



P 



'nq 



(17) 
(18) 



It is obvious that estimates (16)-(18) are non-trivial only when the expression in the 
brackets is positive. Let p < 1/2, nq — > 00 and q — >• 0. Combining Theorems 3.1 and 3.6, 
for sufficiently large n, we obtain 

C 4 q<\\F n -HD*\\<C 5 q, 

cJ-< \\F n - HDlWcc <C 5 J^-, 
V n V n 

< \Wn - HD^\\ W < C b q^. 

Of course, the last estimate, as well as the one in (18), is of interest only if qy/n!} — > 0. 
Similar results can be obtained for the signed CP approximations. 

Theorem 3.7. Let condition (8) be satisfied and let nq>l. Then, for some absolute 
constants Cq and Cj, 



\\F n - H exp{xx(G - I)}D%\\ > cJ £ ( \%\ - C 8 ^±J 

\F n - ffexp{xi(G - /)}^ l ||oo > — ( |7s| - C 8 ^±J 

n V v n 9 



li^-^expixxCG-J^Hw^Csg pys»|-C, 



(19) 
(20) 
(21) 



Let ng — > 00 as (7 — > and p — > p. Also, assume that n is sufficiently large so that the 
right-hand estimates of (19)-(21) are positive. We then have 

C 8 W| < 1 1 F n - H exp{x! (G - 1) }D% 1 1 < C 9 Jl, 

— ^WFn-HcMMG-mDnoo^—, 

n n 



C 8 q < \\F n - 7Jexp{xi(G - I)}D n 2 \\ w < C g q. 



3.5. Asymptotically sharp constants 



In the previous section, we proved that upper and lower bound estimates are of the same 
order, provided that nq is large and q is small. As it turns out, if, in addition, p is small, 
then it is possible to obtain asymptotically sharp constants. 
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Theorem 3.8. Let p < 1/4, q < 1/30 and nq > 1 . Then 



where 



\F n ~HDi\\-A n \<Cq[p + q 



/nq 



\F n -HDt\\ 00 -A 12 \<C x /?-(p + q+-^= 

n \ ynq 



1 



\F n -HDi\\ w -A 13 \<Cqy/riq [p + q 



An = 



4| 72 | 



1 72 1 



-4 



13 



7igV27te 7i yHiy/Znnq q^/lYM 

As a consequence of (22), we note that if p — > 0, q — > and — > oo, then 

||F n -^||- ^ 



'27te 



(22) 

(23) 
(24) 



Similar relations can be obtained for the local and Wasserstcin norms as well. 



4. Applications of Markov binomial models 

In this section, we discuss some areas where the results of our paper can be applied. 

(i) Aggregate claim distribution in the individual model. Consider a portfolio of n 
risks. Each risk produces a positive claim amount during a certain reference period. The 
aggregate claim of the portfolio is then 

S hld = X 1 +X 2 + --- + X n . 

It is usually assumed that all Xj are independent. However, the independence of claims 
does not always reflect reality. For example, an accident involving a tourist group, life 
insurance for a husband and wife or pensions for workers of the same company are likely to 
produce dependent risks. For discussion of the dependence of risks and further examples, 
see [10, 16] and [26]. 

Compound Poisson and signed compound Poisson approximations in the independent 
case of an individual model have been quite thoroughly investigated; see, for example, 
[15, 17]. On the other hand, there are only a few results for the total variation metric 
for dependent risks. Dhacne and Goovaerts [10] investigated a similar model (although 
not explicitly Markovian) under an assumption which, in our notation, is equivalent to 
q m = 0. However, under such an assumption, one cannot expect the limiting law to be 
compound Poisson. Therefore, we have excluded this peripheral case from this paper, 
assuming q to be small, but not identically zero. In [16], Poisson approximation in the 
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general setting of dependent risks was discussed. However, in our case, their result is not 
applicable since for small q, the distribution of the approximated sum is not close to the 
Poisson distribution, but rather to the compound Poisson law. 

Let us assume that aggregated claim amount S md of the portfolio consists of N inde- 
pendent groups of risks. We assume a homogeneous model for each group of risks with 
Markovian dependence. Let each risk have a two-point distribution. More precisely, let 

N n m 
m— 1 j — 1 

Here, X™ and X l h are independent if m ^ I. We assume that each risk of the mth group 
can produce a claim of size a m . Moreover, the dependence of risks of the same group is 
Markovian: P(Xf = a m ) = q m , P(X[ n = 0)=p m and 

Ppff = a m \X™_ x = 0) = q m , P(Xj" = 01XJL, = 0) = p m , 
P(X™ = a m \X™ 1 = a m ) = p m < 1/2, P(Xf = 0\X™ x = a m ) = q m , 

Pm + q m = c[ m +P m = 1 i Pm,5 m G(0,l), m = l,2,...,N,j = 2,...,n m . 

The results of the previous sections can now easily be applied. We illustrate this with 
just one example. Let us define a compound Poisson variable in the following way: 

N N m 

gcp _ ^ Yj m . 

m— 1 j— 

Here, Yj m are i.i.d. geometric random variables, P(Yj m = k) = qmP^T 1 , k — 1, 2, . . . , and 
N rn is a Poisson random variable with parameter n m q m q m /(q m + q m ). The random vari- 
ables N m , m = 1, 2, . . . , N, are independent and also do not depend on Yj m . Denote the 
distributions of S md and S cp by F lnd and F cp , respectively The characteristic function 
of F cp is then given by 



F cp (t)=expi 



n m q m q m (e ltam - 1) 



(?m +q m )( 1 ~p m e ita ™) 
Also, we have the following estimate of approximation: 

N 



l^iBd _ F c P || < c J2 [q m {jp m + q m )mm{\, {n m q m )-^) 



m—1 

+ min (9m' n mlm) + (Pm + <7 m )e" Cin " 



(25) 



Note that the approximation is closer if all q m are small. 

For the proof of (25), one should use the triangle inequality, thus reducing the problem 
to N estimates of Markov binomial distributions concentrated on 0, a m , 2a m , . . . . The 
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total variation metric is invariant with respect to norming. Therefore, without loss of 
generality, one can switch to integer numbers and apply (9) N times with po = 0. 

It is obvious that the second order estimates and estimates in Wasscrstein metric can 
be obtained in a similar way. 

(ii) System failure models. The Markov binomial distribution naturally arises in 
weather and stock market trends. It is also a natural model for system failure situations. 
As an example, we present one model from Sahinoglu [29], who considered an electric 
power supply system with operating and non-operating states throughout a year-long 
period of operation, discretized in hours. Let Mi be the margin values at hourly steps, 
that is, 



where TPG denotes total power generation, Li denotes power demand (hourly peak load 
forecast) and X denotes unplanned forced outages. Let Yi be an indicator of {Mi < 0}. 
Then S = Y\ + I2 + ■ ■ ■ + Y n represents cumulated hours of negative-margin hours, that is, 
the unavailability of power at the nth hour. It is natural to assume that S has a Markov 
binomial distribution. Notably, the Markov chain Yi, Y2, • ■ • is non-stationary. Therefore, 
many known results about the compound Poisson approximation cannot be applied di- 
rectly. Also, the results of our paper relax the assumptions on transition probabilities 
from Sahinoglu's model and give estimates of the accuracy of approximations. Further, 
as shown in Sahinoglu [29], page 49, the probabilities of the compounding geometric law 
under certain assumptions can be viewed as probabilities for the number of trials required 
to repair the system. 

(hi) Industrial applications: sampling plans. A basic assumption in standard accep- 
tance plans for attributes is that the characteristics of items in the lots are i.i.d. Bernoulli 
variables. Recently, however, the focus has been on monitoring the ongoing production 
process by inspecting the items sequentially. In such cases, the quality levels of suc- 
cessive items arc statistically dependent and it has been found in practice that the 
Markov-dependent model is a very useful one; see [24]. Indeed, Bhat et al. [3] modi- 
fied the standard acceptance sampling plans and proposed sequential single sampling 
plans for monitoring Markov- dependent production processes. Vellaisamy and Sankar 
[35] proposed optimal systematic sampling plans for Markov-dependent processes. We 
will outline some possible new research directions in this field. 

5. Auxiliary results 

We now introduce further notation: 



Mi = TPG -X-L. 



0-1=11, 



a-2=l2 + — , 



«3 = 73 + aia 2 



3 
1 



(26) 



3 



OO 



OO 



Y = G 



I, B = J2(ph 



qiy, K = Y,(jPh-qI-2 lx Y)\ (27) 



j=0 
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In the following two lemmas, C(k) denotes an absolute positive constant depending on 
k. Throughout this paper, we set 0° = 1. 

Lemma 5.1. Let t > 0, k € {0, 1, . . .} and < p < 1. Also, let M be a finite (signed) 
measure concentrated at Z. Then, for Y defined in (27), 



\Y 2 e tY \\ < 



tc 



\Y k e tY \\ < 



(i) 



fc/2 



(29) 



If p< 1/2, the 



|y fc e ty iu< c{k) 



t (fc+l)/2' 



|yM||w>|||M||, \\YM\\>^\\{h-I)M\\. (30) 



Proof. The estimates in (29) follow from the properties of the total variation norm and 
results in [27] and [11]. The first estimate in (30) is a consequence of the inversion formula 
and the following inequalities: 



ReY(t) < 



sin - , 

1+p 2' 



in*)i< 



sin ■ 



Here, Re{-} means the real part of the complex number. In view of the relation between 
total variation and Wasserstein norms (see the Introduction), we get 



||MF|| W = 



(i 1 -i)mJ2p'Ij 



w 



3=0 



\\M{I 1 -I)\\ = \\MY(I-pI 1 )\\<\\MY\\{l+p), 



I Ml 



3=0 



< 



3=0 



The results in (30) now follow easily. 

For our asymptotically sharp results, we need the following lemma. Set 
1 d fc 



□ 



fk(x) 



2ndx k 
(ft = 0, !,...)• 



-x*/2 



<Pk\\i 



tp k (x)\dx, H^fclloo = sup |cp fc (a;)| 
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Lemma 5.2. Let t > and k = 0, 1, 2, PPe i/ien /iawe 

\\(h - I) k e^- 



IMIi 


< 


C(k) 




t (k+l)/2 


ll^fe oo 


< 


C{k) 


^(/c+l)/2 


tk/2+1 ' 


l!vfc-i||i 


< 


C(k) 


t (fe-l)/2 


t k/2 



(fc^O). 

The proof follows from a more general Proposition 4 in [28] . 
Lemma 5.3. If N > and < a <p < 1, then (I + aY) N is a CP distribution. 
Proof. Note that 

(I + aY) N = cxp{-7Vln(l - a)(F - I)}. 
Here, F is a distribution concentrated on {1, 2,3,.. .} with 

The last relation obviously completes the proof. 



□ 



Before we proceed to our main lemma, we need some additional facts about F n . Simi- 
larly to [5] (see also [7]), it is possible to check that under assumption (8), we have 



where 



Ai, 2 (t) 



Wi, 2 (t) 



F n (t) = A?(t)Wi(*) + A5(t)W a (t), 
pe [t +p±D 1 / 2 {t) 



(31) 



Po 



1± 



q + q+pi&t 1-po ( q + q + (2q-p)(e it -1) 



£>V3(t) 

D(t) = (pe H +p) 2 +4e H (q-p). 
This allows us to write F n as 

F n = A?Wi + A!} Wa 
and to express Ai :2 and Wi,2 as the following series: 



1± 



DV2(i) 



A 1 = / + ai r + i{(i + g)/-p/ 1 + 2a 1 y}^ 1 / 2 ^(-i)^^, 



(32) 



(33) 



1128 V. Cekanavicius and P. Vellaisamy 

A 2 =p7 1 -g7+(7-A 1 ), (34) 

Wi, 2 = \ |i ± [(<? + 1) i + p(/i - /)]*: £ ( " J /2 ) (-i) J : i j | 

(35) 

±{l-p ){q-p){h-I)KY j (- l ' 2 \{-iyV. 



The following lemma is used as the main tool in the proofs. 



Lemma 5.4. 7/ condition (8) is satisfied, then 

3 

Ai =I + ^fl/ J ' +C*g 3 (p + g)y 4 9, (36) 
3=1 

3 

]nk 1 =^~ 1j Y i +Cf{p + q)Y i ®, (37) 

3=1 

19 

inA 1 = 7l r + — 7l y 2 e, (38) 

60 

l|A 2 ||<^, HAx-711^0.1, (39) 

Wi =7 + (xi + ^ 2 )y + Cg(p + 9)r 2 e, (40) 

Wa = C(p + g)(Ii --f)6, ||W 2 ||<7. (41) 
Tor any finite signed measure M on Z and any £ > 0, u;e Ziawe 

||Mexp{ilnAi}|| < C||Mexp{(i 7l /30)y}||, (42) 
||M7Jj|| <C||Mexp{(t 7l /30)y}||, j = l,2,3. (43) 

Estimates (42)-(43) also hold for the local norm. 

Proof. We have 

«i=7i<^, -^—<j^—<2, ||r||<j|G||+l = 2, 
30 q + q 1 - p 



\\K\\<J2(p + q + 4 ai y <3, (44) 

3=0 



||7|| <9-4-r -4( l + -^-Jl+p) ) <0.4. (45) 
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Note that 



We have 



1/2 
2 



E 

3=1 



1/2 
3 



1 

16' 



1/2 
3 



- 128' J ~ 



1/2 
j 



< I + M + M! + AM! < 0.5642 
11 - 2 8 16 128 0.6 - 



and 



[7(1 + g) - p/x + 2 ai F]L = 4 * [g 2 / + p(g + «)(/- p/i)]if 



jfy 2 

"(<Z + ?) 2 



Consequently, 



Ai = i + 7i y + io.5642 7 iy 2 e = / + 1.282171^0 = /+ o.ie 



and 



~ (-1V+ 1 
In Ax = Ai - / + J2 ~ : ( A i - T Y 



3=2 



71 



1 00 19 

y + o.282i 7 iy 2 e + -i.282i 7 2 y 2 ^(o.i)^ 2 e = 7i y + —71 y 2 ©. 



J=2 



Moreover, ||A 2 || <p+g + ||Ai - I\\ < 19/30. Thus, we have proven (38). We use this 
estimate for obtaining (42). By the properties of the total variation norm, we have 



I Me 



tin A, 



< 



Mexp<; — y 

^ 30 



cxp 



[29t7i y 19*71 y2Q 



30 



60 



Applying Lemma 5.1, we prove that the second norm is majorized by 

^ 1 19 , f 29 1 '' ^ e- /57r\ r ^ 

l+> — — t7iy 2 exp<^ *7iY> <l+> ;= — <c. 

60 ' M30r n J " ^ ^v^nr V 58e J " 

r— 1 v y r=l v v ' 



The last two estimates obviously lead to (42). The estimate (43) is proved similarly. For 
the proof of (36), note that 

Ai = I + iiY - ~^--Y 2 b 2/ + P(9 + Ph)]K + CtY A S, (46) 

(q + q) 2 
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h-I = Y(I- ph), (q + q)B = I + p(h - I)B, 
h-I= 2YQ, (q + q)(I-ph)B = ql-pq(l - ph)BY, 



B 



1 



pq 



p 2 q 



;Y(h-I)B, 



q + q (q + q) 2 (q + g) 2 
K = B + 2jiYKB = B- 2^ X B 2 Y + A^\B 2 Y 2 B 2 K. 



(47) 
(48) 

(49) 

(50) 



Substituting (48)-(50) into (46), we obtain (36). Taking into account (39), we obtain 



In 



^ f-1 V' +1 

Ai = V M — (Ai - iy + c(Ai - /) 4 e. 

H 3 



Now, for the proof of (37), it suffices to use (36). From (50) and the first relation in (48), 
we get p{I\ - T)K=p{h - I)B + CpqY 2 Q. Moreover, 



p(h - I)B 



^ zY _^h-t)BY = P3L + C pqY 2 e. 



q + q 



q + q 



q + q 



The last two equations and (35) allow us to prove (40). Since W\ + W2 = I, we easily 
obtain the first relation in (41). Now, 



\\wh< 



Ul + (q + q)3jr 



-1/2 
3 



OA 1 



\q-p\-2-3j2 

3=0 



-1/2 
J 



OA 1 < 7. 



Thus, Lemma 5.4 is proved. For the lower bound estimates, we need the following result. □ 
Lemma 5.5. Let M be concentrated on Z, aeM and b > 1. Then, 



\\M\\>C 
|Af|U>? 



-ita 



dt 



' t2 ' 2 M 



-ita. 



dt 



(51) 
(52) 



The estimates (51) and (52) remain valid if e 1 I 2 is replaced by te * / 2 



Lemma 5.5, with |jAf|| replaced by the uniform norm of M, was proven in [32]. Since 
the uniform norm is majorized by the total variation norm, (51) also holds. 



Lemma 5.6. If (8) is satisfied, then, for all \t\ < tc, 

I exp{n^i(r(i) - it/q)} - 1| < Cnqt 2 , 
\D%{t)\ < 1, |i3^exp{-iin7i/g}- 1| < Cnqt 2 . 



(53) 
(54) 
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Proof of Lemma 5.6 is straightforward and therefore omitted. 
Finally, let us introduce an inverse compound measure for H . Let 




Lemma 5.7. If (8) is satisfied, then 



(55) 



and, for any (signed) finite measure M concentrated at Z, 



\\MH\\ >c~ 2 ||M| 



||Mexp{-p o7l y}|| > \\M\ 



(56) 



The estimates in (56) remain valid if the total variation norm is replaced by the local 
norm. 

Proof. Estimate (55) easily follows from the property ||e|| M < e" M " ; see the Introduction. 
Now, = WMHH^W < WMHWWH-m < c 2 \\MH\\. Since cxp{p 7iY"} is a distri- 

bution, its total variation is 1. Therefore, ||M|| = ||Mexp{— pojiY} exp{po7i^}|| < 
||Mexp{— Po7i^}|| ■ Estimates for the local norm are proved similarly. □ 



For upper bound estimates, we use an adaptation of Le Cam's [23] approach which deals 
with convolutions of measures. 

Proof of Theorem 3.1. Without loss of generality, we can assume that (8) holds. We 

have 



6. Proofs 



\\F~n - HD*\\ < ||A» - D^WWW^l + \\D^(W 1 - H cxp{~ Poll Y})\\ + ||A 2 ||"||PF 2 ||. 



Further, in view of Lemma 5.4, 




D?\\ < D 



/ (exp{T[rtlnAi - 717^]})^ dr 



Jo 



< n 




|| [In Ai -7iF]exp{TnlnAi + (1 - T)n7iF}|| dr 



< Cn 1 1 [In Ai - 71 Y] exp{ (n 7 i /30) Y} \ \ < Cnq 2 \ \ Y 2 exp{ (n 7i /30)Y } \ \ . 
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By Lemma 5.4, 

Wi - Hexp{-p oll Y} = [W X -I- (xi + x 2 )Y] + [I + (xi + x 2 )Y - (I-p o7l Y)H] 

+ [H(I - Poll Y - cxp{- Po7l y})] = Cq(p + q)YQ. 

Taking into account the last two estimates, applying Lemma 5.1 and estimating HW2II 
and HA2II by (41) and (39), we complete the proof of (9). The estimates in (10) and (11) 
are proved similarly. □ 

Proof of Corollary 3.2. Following the proof of (42), one can prove the same property 
for Hi. Also, 

\\HDi — HHi\\ < CX\\(Di — Hi) exp{(n 7 i/30)y}|| < Cnq 2 \ \ Y 2 exp{ (n 7 i /30) Y } 1 1 . 
The rest of the proof is obvious. □ 
Proof of Theorem 3.2. We have 

\\F n -HD^{I + n l2 Y 2 )\\ 
<||A 2 ||"||^ 2 || + ||^ 1 ||||Ai l -^!| 

+ 11^11 ||I? r 1 l (c n721 ' 2 -I-n l2 Y 2 )\\ + \\D^{I + n l2 Y 2 ){Wi-Hc~^ Y )\\. 
Similarly to the proof of Theorem 3.1, and using (43), we obtain 



|AJ-DJ|| <Cn 



In 



Ai-7iF-72^ 2 ] / cxp{rnlnAi + (1 - T)[nyxY + n-f 2 Y 2 }}dr 
Jo 



< Cn\\ [In Ax - 7l y - l2 Y 2 ] cxp{(n 7l /30)y}|| 

< Cnq 2 {q+p)\\Y 3 cxp{(nj 1 /30)Y}\\ 



and 



f 1 

(n l2 Y 2 ) 2 D^e Tn ^ Y (l-r)dr 
Jo 



\\D?(c n ^ Y2 -I-n l2 Y 2 )\\ < 

<C(n 72 ) 2 ||r 4 cxp{(n 7l /30)r}||. 

Note that for any signed finite measure M, 

\\D^{I + n l2 Y 2 )M\\ <||^i /2 ||M(l + 7i| 72 |||r 2 ^ /2 ||)<C|| J D' 1 l/2 M||. 

The rest of the proof is very similar to the proof of Theorem 3.1 and is hence omitted. □ 

Proof of Theorems 3.3, 3.4 and 3.5. The proofs are very similar to those of Theo- 
rems 3.1 and 3.2. From Lemma 5.4 and the definition of the exponent measure, it is not 
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difficult to show that 

Wi - e* lY H = [Wi - 1 - (xi + x 2 )Y] + [I+ (xi + x 2 )Y - (I + x x Y)R\ 
+ H(I + xiY - c* lY ) = Cq(p + q)Y 2 0, 
||A?Wi -D^Hc xlY \\ < ||A?-£>J||||Wi|| + \\D 7 2 l (W 1 - Hc" lY )\\. 

Now, it is not difficult to prove Theorem 3.3. Theorem 3.5 is proved similarly. For the 
proof of Theorem 3.4, one should use Theorem 3.5, the triangle inequality and the fact 
that 

!|#»(/ + n 73 r 3 )-^!| = D% / (l-r)e™^ y3 (n 73 y 3 ) 2 dr 

Jo 

<C(n 73 ) 2 ||r 6 exp{n 7l y 2 /30}||. 
For the last estimate, we have used the same argument as in the proof of (43). □ 
Proof of Theorem 3.6. Taking into account Theorem 3.2, (30) and (56), we get 

V 



\F n - HDlW > n\ l2 \\\HDiY 2 \\ - Cql q 



>C w n\ l2 \\\D1(h-iy\\-C lx q [q 



P 



'nq 



(57) 



Let z — t/(hs/nq) and /i = wyx/q. The constant h > 1 will be chosen later. Applying 
Lemma 5.6, we then obtain 



J = 



> 



e~ t2 / 2 z 2 dt 



Here, 



Jt= [ e-*/*2?\D?(z)e- i " t -l\dt<Crq f z^ ' 2 &t ■■ 

JM JTg. 



Jl-J2 



c 



(58) 



h A nq' 



J 2 = / e 



\D^{z)e- lz ^\\(e lz - l) 2 - {iz) 2 \ dt< 



C 



h?nq^fnq 



Combining the last two estimates with (58) and choosing h to be a sufficiently large 
absolute constant, we obtain 



J> 



C12 
h 2 nq 



1 - 



C13 C 



14 



h 2 hy/nq 



> 



C 



15 



nq 



Applying Lemma 5.5 and substituting the result into (57), we get (16). Estimates (17) 
and (18) are proved similarly. 
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For the proof of Theorem 3.7, one should use Theorem 3.4 and take iexp{— 1 2 /2} 
instead of exp{— 1 2 /2}. The proof is then almost identical to that of Theorem 3.6 and is 
hence omitted. □ 



Proof of Theorem 3.8. We have 
\WFn-HDi\l-Au\ < \\F n -HD^(I + n l2 Y 2 )\\ 

+ \\(He-P°^ Y - I)D?nj 2 Y 2 \\ + n\j 2 \ 

, n \l2 



Y 2 --(h-I) 2 )D? 



(h - If D n x - exp 



{>""}) 



- An 



One should now apply Theorem 3.2, (48), Lemmas 5.2, 5.1 and the following, easily 
verifiable, relations: 



' \ " i q q z 



and 



q "oV« 



D 1 -exp\ ^(h-I) 



: exp /^(A - 1) | ^exp|^(Ji - I) 2 e| - = Cpq(h - I) 2 Q. 



Note that 



(h If (^i - exp | ^ (h I) | ) £ D?"' exp 



All other estimates are obtained similarly. 



□ 
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